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Csi ! Abstract 



In this article we derive the integral equation that controls the momentum dependence of the ef- 
fective gluon mass in the Landau gauge. This is accomplished by means of a well-defined separation 



D , of the corresponding "one-loop dressed" Schwinger-Dyson equation into two distinct contributions, 



one associated with the mass and one with the standard kinetic part of the gluon. The entire 



^ ' construction relies on the existence of a longitudinally coupled vertex of nonperturbative origin, 

^O . which enforces gauge invariance in the presence of a dynamical mass. The specific structure of the 

a^ 

C^ I resulting mass equation, supplemented by the additional requirement of a positive-definite gluon 

f^ ■ mass, imposes a rather stringent constraint on the derivative of the gluonic dressing function, which 

is comfortably satisfied by the large-volume lattice data for the gluon propagator, both for SU{2) 

and SU{3). The numerical treatment of the mass equation, under some simplifying assumptions, is 

$H ' presented for the aforementioned gauge groups, giving rise to a gluon mass that is a non-monotonic 

C^ ■ 

function of the momentum. Various theoretical improvements and possible future directions are 
briefly discussed. 
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I. INTRODUCTION 



A large body of recent high-quality lattice results indicate that the gluon propagator 
and the ghost dressing function of pure Yang-Mills theories, computed in the conventional 
Landau gauge, are infrared (IR) finite, both in SU{2) [l|-[5J and in SU{3) |6H9|. These 
important results have sparked a renewed interest in the important issue of dynamical mass 
generation in non-Abelian gauge theories, and especially in QCD lOMlSl]. Specifically, as has 
been suggested in a series of works, the finiteness of these quantities may be interpreted as 
a direct consequence of the generation of a non-perturbative (momentum-dependent) gluon 



mass, which acts as an IR cutoff of the theory ly, [l7| . In the picture put forth in these 
articles, the fundamental Lagrangian of the Yang-Mills theory (or that of QCD) is never 
altered; the generation of the gluon naass takes place dynamically, without violating any of 
the underlying symmetries [l2|. 
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16|. 



Given the non-perturbative nature of the mass generating mechanism, its study in the 
continuum proceeds through the Schwinger-Dyson equations (SDEs) that govern the dynam- 



ics of the various Green's functions of the theory 11 
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42011 . and especially of the gluon 



propagator, A(g^). The main conceptual and technical challenge in this context is to obtain 
as a solution of these integral equations an IR-finite gluon propagator [i.e., A^^(O) = m^(0)], 
without interfering with the gauge invariance (or the BRST symmetry) of the theory, en- 
coded in the Ward identities (Wis) and Slavnov- Taylor identities (STIs) satisfied by the 



Green's functions under study lll-ll3l|. A self-consistent framework for enforcing the crucial 



property of gauge invariance at the level of the truncated SDEs is provided by the synthesis 
of the pinch technique (PT) [jj, [id, UM^ with the background field method (BFM) [24 1. 



In the presence of a dynamically generated mass, the (inverse) Euclidean gluon propagator 
assumes the form A~^(g^) = q^J{q^) + m?{q^), where the first term corresponds to the 
"kinetic term" , or "wave function" contribution, whereas the second is the (positive-definite) 
momentum-dependent mass [l5|]- However, to date, practically all studies attempting to 
determine the IR behavior of the gluon propagator from SDEs eventually boil down to the 



16|, 



solution of some integral equation involving the entire gluon propagator A(g^) 12, llJ, 
rather than its two components, J{q^) and rn?{q^). This is to be contrasted to what happens 
in the analogous studies of chiral symmetry breaking, where one derives a system of two 
coupled equations, one determining the "wave function" ( "kinetic part" ) of the quark self- 



energy, and one determining the dynamical (constituent) quark mass |25l . |26| . Of course, 
in the case of the quark self-energy the above separation of both sides of the corresponding 
SDE (quark gap equation) is realized in a direct way, due to the distinct Dirac properties 
of the two quantities appearing in it, while in the case of the gluon propagator no such 
straightforward separation is possible. However, a well-defined procedure, first outlined 
in [l5|, and explained here in more detail, allows for an analogous separation even in the 
case of the gluon propagator. The purpose of the present article is to identify and isolate from 
the SDE of the (Landau gauge) gluon propagator the dynamical equation that determines 
the evolution of the gluon mass, study its main properties, and find approximate solutions 
for m?{q^). 

As has been emphasized in some of the literature cited above, a crucial condition for the 
realization of the gluon mass generation scenario is the existence of a longitudinally coupled 
vertex, to be denoted by V , which must be added to the conventional (fully-dressed) three- 

[ ! 

gluon vertex, denoted by F [l5i] . Specifically, the vertex W = T + V satisfies the same STIs 
as r, but now replacing the gluon propagators appearing on their rhs by a massive ones 
(schematically, A — )■ A^)- The dynamical reason for the emergence of this special vertex, as 
well as its diagrammatic realization in terms of Feynman graphs, is intimately connected to 



the well-known Schwinger mechanism 27|,|28|, which enables the non-perturbative generation 



of a gauge-boson mass. In particular, one assumes that the strong QCD dynamics give rise 



to longitudinally-coupled composite (bound-state) massless poles |29l-l35|. These poles act 
like Nambu-Goldstone excitations, in the sense that they preserve the form of the STIs of 
the theory in the presence of a mass, but they are not associated with the breaking of any 
local or global symmetry. 

It turns out that the way the vertex V generates the mass at the level of the SDE is 
by introducing a "deviation" from the so-called "seagull identity" [given in Eq. (14. 3p ]. The 
role of this identity is to enforce the masslessness of a gauge boson (gluon or photon) when 
massive propagators appear inside its loops, assuming always that the WI and STI's are 
maintained, i.e., the transversality of the (gluon or photon) self-energy is preserved. For 



example, as explained in 15|, in scalar QED it is exactly this identity that enforces the 
masslessness of the photon at the level of the "one-loop dressed" SDE; in this case the 
massive propagator entering into the loop is that of the charged scalar field. The crucial 
point is that if the "massive" STI were to be enforced by only modifying F [i.e., by carrying 



out the replacement A — )• Am in the closed expressions obtained for T by solving the STIs 
it satisfies, see Eq. (IA3|) ]. then the "seagull identity" would force the (would-be) gluon mass 
to vanish, i.e., would lead to the invalidation of the entire mass generation mechanism. 
The fact that the missing part for satisfying the "massive" STI is instead provided by the 
longitudinally coupled V has the far-reaching consequence of finally furnishing a non-trivial 
equation for the mass. Thus, the equation for the gluon mass is determined as the amount 
by which the seagull cancellation is distorted due to the presence of the vertex V. 

To be sure, one could in principle determine the closed form of V by resorting to a 



procedure similar to that described in 36|] for F, namely write down the most general 
structure allowed for a longitudinal vertex with three Lorentz indices, and then determine 
the corresponding form factors from the WI and STI that this vertex must satisfy. It turns 
out, however, that V enters into the SDE for the Landau gauge gluon propagator (in the 
PT-BFM scheme) in a very particular way, which renders its closed form unnecessary; all 
one needs for the derivation of the mass equation is to postulate the existence of V {i.e., 
assume that it is not identically zero) and that it satisfies the required WI and STIs. 

In principle, the mass equation obtained in Eq. (I5.23P must be accompanied by the 
corresponding equation determining the kinetic term J{q'^); the solution of the resulting 
system of two coupled integral equations will then furnish the behavior of m?{q^) and Jiq^), 
and therefore that of A(g^). The technical limitation in realizing these procedure is the 
dependence of the equation for J{q^) on the various form factors comprising the ghost-gluon 
kernel; the latter enters into play through the form of the vertex W [see Eqs. ( 1A3I) -( IA4|) ] . 
The way to circumvent this problem is to actually solve Eq. (I5.23p for m?{q^) using as input 
for the A appearing in it the available lattice data, both for SU{2) [1] and SU{3) |6|. 

It turns out that the specific form of the mass equation in Eq. (I5.23P introduces a 
non-trivial constraint on the precise behavior that A must display in the region between 
(1-5) GeV^. Specifically, in order for the gluon mass to be positive definite, the first deriva- 
tive of the quantity g^A(g^) (the "gluon dressing function") must furnish a sufficiently 
negative contribution in the aforementioned range of momenta. Interestingly enough, the 
A obtained from the lattice has indeed this particular property. This is to be contrasted 
to what happens, for example, in the case of a simple massive propagator l/(g^ -|- m"^) or 



with the Gribov-Zwanziger propagator q^ /{q^ + rn^) (with m constant) 37 



the deriva- 



tives of the corresponding dressing functions, q^ /{q^ + rn^) and q^jiff + m'^\ respectively. 



are positive in the entire range of (Euclidean) momenta, thus excluding the possibility of a 
positive-definite gluon mass. 

The article is organized as follows. In Section |Tl] we introduce the necessary notation and 
review briefly the aspects of the PT-BFM formalism relevant to this work. In Section IIIII 
we explain in detail the modifications that must be introduced to the three-gluon vertex of 
the theory in order to treat the generation of a gluon mass in a gauge invariant way {i.e., 
preserving the STIs of the theory). In particular, the importance of the nonperturbative 
vertex V and its special properties are emphasized, and the changes introduced to the 
corresponding SDE during the transition from massless to massive solutions are discussed 
in detail. In Section |IV] we outline the precise criteria that will lead to the separation of 
the SDE for the gluon propagator into two equations, one for the kinetic part and one 
for the mass. The central role of the "seagull-identity" in carrying out this separation is 
stressed, and some explicit characteristic calculations are presented. Then, in Section |V] we 
combine the ingredients introduced in the previous sections and derive the final form of the 
dynamical equation for the gluon mass in the Landau gauge. In Section |VT] we first study the 
implications of the gluon mass equation in the limit of vanishing physical momentum. Then 
we solve an approximate form of this equation, using lattice data as input for the "unknown" 
quantity A. The solution for the gluon mass so obtained is then appropriately "subtracted 
out" from A, giving rise to an estimate for the quantity J{q^)- These ingredients are then 
combined to construct the renormalization-group (RG) invariant gluon mass, appearing 
in the usual definition of the effective QCD charge within the PT-BFM framework. Our 
conclusions and discussion of the results appear in Section IVIII Finally, some technical 
points are presented in the Appendix. 

II. GENERAL FRAMEWORK 

In this section, we set up the necessary notation and review some of the most salient 
features of the PT-BFM framework, putting particular emphasis on the form of the SDE 
for the gluon propagator, and the various field-theoretic ingredients appearing in it. 

The (full) gluon propagator A^^(g) = 6"'^A^^{q) in the renormalizable R^ gauges is defined 
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FIG. 1: The "one-loop dressed" gluon contribution to the PT-BFM gluon self-energy. White 
(respectively, black) blobs represent connected (respectively, 1-PI) Green's functions; a gray circle 
on the external legs indicates background gluons. Notice that within the PT-BFM framework these 
two diagrams alone constitute a transverse subset of the full gluon SDE. 
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with 



p,MW) + ^ 
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A;J(g) = ^ [P,M)^-\q') + {VOwA , (2.1) 



P^Luiq) = Qfiv - -^ 



(2.2) 



the diinensionless transverse projector, and ^ the gauge fixing parameter. The scalar cofactor 
A(g^) appearing above is related to the all-order gluon self-energy Ii^u{q) = P^^{q)Il{q'^) 
through 

A-\q')=q' + tU{q'). (2.3) 

In addition, it is convenient to define the dimensionless function J(g^) as [39| 

A-\q') = q'j{q'). (2.4) 

Evidently, J(g^) coincides with the inverse of the gluon dressing function, frequently con- 
sidered in the literature. 

The starting point of our dynamical analysis is the SDE governing the gluon propagator. 



Within the PT-BFM framework that we employ [lOHlTl. |21|-|24|. one may safely truncate 
the SDE series down to its "one-loop dressed version" containing gluonic contributions only, 
given by the diagrams (ai) and (02) shown in Fig. [1] Specifically, due to the special Feyn- 
man rules of the PT-BFM, and in particular the QED-like Ward identities satisfied by 
the fully-dressed vertices, gauge invariance remains exact, in the sense that the resulting 



(approximate) gluon self-energy n^,y(g) is still transverse, i.e., 

q-'U^M = 0. 



(2.5) 
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FIG. 2: The BQQ three-gluon vertex. 



The PT-BFM equation for the conventional propagator reads, in this case, 



[1 + G{q')] 



2^12 



where 



1 
2 



l«iV = 7^9'Ca lf'^;;Uq,k,-k - q)A"^{k)A^''{k + q)r,,^{q,k,-k - q) 



(«2) 



fiu 



g'c. 



g,jA';{k) + (1/^-1) JA,M 

Jk Jk 



(2.6) 



(2.7) 



with Ca being the Casimir eigenvalue of the adjoint representation [Ca = N for SU{N)], 
and the d- dimensional integral measure (in dimensional regularization) is defined according 
to 

The vertex F is the fully dressed version of the trilinear vertex involving one background and 
two quantum gluons {BQQ vertex for short, see Fig. [2]); at tree-level (all momenta entering) 



rLt(?'^'P) = rLt(?'^'P) + (i/OC(?,^,p) 



(2.9) 



with 



rS.(g,r,p) 



9f,uir - p)a + gau{p - q)^l + QcxM - r)u, 

9t,v{r - p)c + gav{p -q + r/Cj^ + 9aM -r- p/Oi 



(2.10) 



Finally, the function G{q'^) appearing in (12. 6 p is of central importance in this entire 
formalism. It is defined as the scalar co-factor of the g^^^, component of the special two-point 
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H^fXq,P,r) = 5,,^ + 




FIG. 3: Diagrammatic representation of the functions A, H and, for later convenience, H. Gray 
blobs represent 1-PI kernels with respect to vertical cuts. 

function A^^{q), defined as (see also Fig. |3]) 

KM = -ig^CA fA';^{k)D{q-k)H,^{-q,q-k,k) 
Jk 



9,Mq') + ^-^Hq 



(2.11) 



where we have introduced the ghost propagator D"-'^{q^) = 6°'^D{q'^), which is related to the 
ghost dressing function F(g^) through 



D{q 



2^ F{q') 



(2.12) 



Notice that in the Landau gauge, an important exact (all-order) relation exists, linking 
G(g^) and L^q"^) to the ghost dressing function F^q"^), namely |40l-l43| 



F-\q') = l + G{q') + L{q'). 



(2.13) 



In addition, the function G{q^) participates in a set of BRST-driven identities, known 
as Background-Quantum identities (BQIs) 4^, |45|], obtained within the Batalin-Vilkovisky 



formalism 



46 



47j . These powerful identities relate among each other the three types of gluon 



propagators that appear naturally in the BFM formalism, namely: (z) the conventional gluon 
propagator (two quantum gluons entering, QQ), denoted by A(g^); (m) the background gluon 
propagator (two background gluons entering, BB), denoted by A(g^); and {Hi) the mixed 
background-quantum gluon propagator (one background and one quantum gluons entering, 
BQ), denoted by A(g^). The corresponding BQIs are 



A(g^) = [l + G(g^)]^A(g^) 



A(g2) = [l + G(g2)]A(g2). (2.14) 

Notice that it is the first of these identities that allows the rewriting of the conventional 



SDE into the PT-BFM form ([21]) [ll|- 



For the rest of the article we will study the gluon SDE in the Landau gauge, ^ = 0. The 



limit of Eq. (12. 6 p as ^ — )■ is rather subtle, and has been presented in 12|. The final answer 
is 



/C^^Ar(g), (2.15) 



i=l 



with 



Ar(g) = JP'-'ik) ' ^^'^^^^f A(fc), 



The vertex F appearing above is the fully-dressed PT-BFM vertex studied in detail in 
and which is related to the full BQQ vertex F appearing in the BFM through 

]faMi^(g,r,p) = fo,^,u{q,r,p) + (l/Or^^^(g,r,p). (2.17) 

Evidently, Ta^y{q,r,p) and Fa^,^(g, r, p) differ only at tree level; specifically, one sees imme- 
diately that 

rlt(?,^,p) = rS.(?,^,p)- (2.18) 

In the rest of this paper, we will refer indifferently to both F and F as the BQQ vertex; in 
addition, in order to simplify the notation, we will drop the "tilde" superscript. 

The vertex F satisfies a (ghost-free) WI when contracted with the momentum q^ of the 
background gluon, whereas it satisfies a STI when contracted with the momentum of the 



quantum gluons (r^ or p^). In particular, 

r^r,^,(g,r,p) = ^(r^) \q'j{q^)P^{q)H^,{q,r,p) - p-'Jip^)P!^ip)H^^ip,r,q) 
p''T^^,{q,r,p) = F{p^)\r'j{r^)P;{r)H,^{r,p,q)-q'j{q^)PMHu^iq,P,r) ,(2.19) 



and the function J is related to the conventional one defined in (12.41) precisely through the 
second equation in (I2.14p . namely 



J{q')=[l + G{q')]j{q'). 



(2.20) 



In addition, as shown in Fig. |3l the auxiliary ghost function H is the same as H after 
converting the external gluon leg into a background leg. An explicit form in terms of J, 
J, H and H of the (longitudinal) form factors characterizing this vertex has been obtained 



m 



36| and reported in Appendix |X1 

One may finally use Eq. f l2.4p to re-express the relations (12.191) in terms of the (inverse) 



scalar functions A, i.e., 



g-r„^,(g,r,p) = A-\p')P,,{p) - /\-\r')P,,{r 



(2.21) 



with analogous expressions holding for the remaining two STls of (I2.19p . At this level this 
appears as a simple rewriting, but this form of writing will facilitate the clarification of 
certain conceptual issues that become relevant when dynamical mass generation is turned 
on (see next section). 

III. VERTICES IN THE PRESENCE OF A DYNAMICAL MASS 



In order to generate a dynamical mass without interfering with gauge invariance and the 
BRST symmetry, one must resort to the Schwinger mechanism 27|, l28|] . The general idea 
is to assume that a longitudinally coupled bound-state pole has been formed dynamically, 
which will modify the structure of the full vertices of the theory 29l-l35l|. This modification, 
in turn, will be responsible for obtaining massive type of solutions from the SDE of the gluon 
where this new vertices will be inserted 12|. It is important to be very precise regarding 
the nature and role of the various ingredients that enter in the ensuing analysis. We will 
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FIG. 4: The IT' three-gluon vertex. Thick (red onhne) internal gluon lines indicates massive 
propagators A^, as explained in the text. 

therefore devote this section to the development and elaboration of the various key concepts 
needed. 

From the kinematic point of view we will describe the transition from a massless to a 
massive gluon propagator by carrying out the replacement (in Minkowski space) 



A 



-l^„2\ 



gV(g^ 



Kn\ql = q'JM)-m\q') 



(3.1) 



The symbol Jm indicates that effectively one has now a mass inside the corresponding 
expressions: for example, whereas perturbatively J{q^) ~ Ing^, after dynamical gluon mass 
generation has taken place, one has Jmi.^f') ~ ln(g^ + rn?). As a consequence, since J^ will 
be the main component in the definition of the QCD effective charge ISj, the presence of 
;he mass term in the argument of its logarithm will tame the perturbative Landau pole 
431 . Of course, as q^ — )■ 0, q^Jmiff) — > 0; therefore, if we are to ensure that this 



16 



41 



procedure will give rise to a non vanishing IR value for the gluon propagator, i.e., A~^ (0) 7^ 0, 
we must have that m^(0) 7^ 0. 

From the dynamical point of view, it is clear that the full three-gluon vertex must be also 
appropriately modified |29l-l31|. Specifically, we will consider a new vertex, to be denoted 
by r', and carry out the replacement 



r 



r' = r^ + V. 



(3.2) 



The new vertex F^ is given by the same (fully dressed) graphs that make up the SDE of the 
BQQ vertex F (Fig. Hj); however, now all internal (virtual) fully-dressed gluon propagators 
are massive z.e., in the non-pole part of the vertex SDE we have A — )■ A^. In addition (or 



11 



as a result thereof), F^ satisfies exactly the set of WI and STIs given in (12.190 . but with 
the replacement J — t- Jm throughout. So, the WI becomes 

g„rr(g,r,p) = p\Up')P^^{p) - r'JUr')P'''{r), (3.3) 

and exactly analogous expressions for the remaining STIs satisfied when F^ is contracted 
by either r or p. Note that all other Green's functions, such as H and H, must be replaced 
by the corresponding Hm and Hm, in the same sense as before (but we will suppress their 
'm' subindex throughout); thus, the diagrams defining these two ghost functions, shown in 
Fig. [3l will now contain massive internal gluon propagators. 

On the other hand, the vertex V represents the pole part of F'; it is totally longitudinally 
coupled, i.e., it vanishes identically when contracted by the three transverse projectors 

P"'"(g)P^'^(r)P'^''^(p)K^,(g,r,p) = 0, (3.4) 

and must satisfy the WI and STI of (12.191) . with the replacement k'^J{k) — >■ —rn?{k), e.g., 

q''V^^Aq,r,p) = -m\p')P^,{p) + m\r^)P^,{r). (3.5) 

Exactly analogous expressions will hold for the STIs satisfied when contracting with the 
momenta r or p. 

An explicit example of such a vertex (which, however, we will not use here), has been 



given in 48|, namely 



V^^^{q,r,p) = ^^i:!li^^pP'^(p)rn\p')-^[m\r')-m\q')]P^{q)P^^{r) 
zq r p 

Zr p q 

+ ^^^^^^j^^P''^{r)m\r')-^-^[m\q')-m\p')]P^{p)PP'^iq). (3.6) 
Zq p r 

The totally longitudinal nature of this vertex is manifest^. 

At this point it is clear that the full vertex F' will satisfy the same WI and STIs f l2.19p 
satisfied by the F vertex before the introduction of any masses, but now with the replacement 



^ Note that this vertex is totally Bose symmetric, satisfying p.Sp with respect to all its legs; instead, the 
vertex considered here satisfies an STI with respect to the quantum legs (r and p). 
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A -^ A,„. Therefore, using Eqs. ([32D, (O, and ([33]), one gets for T' the WI 

= \p'Jm{p') - m\p^)]P^^^{p) - [r\Ur^) - m\r^)]P^^^ {r) 

= l^-^\p')P'''{p)-K!{r')P"'{r). (3.7) 

Similarly 

r^r"^^(g,r,p) = F{r^)[^-^\q^)Pi:;{q)H,M.r^P) - K!{P^)P!:{p)H,c.{p.r,q)\, (3.8) 
p.r"^'^(g,r,p) = F{p^) [A-i(r2)P;(r)i^.,(r,p,g) - A„i(g2)P,^(g)/f,^(g,p,r)] . (3.9) 

It is very important to emphasize that, even though the new (massive) WI is obtained 
from the old (massless) one through the replacement A — )■ A^, the new vertex F' is not 
obtained from the old one, F, by means of the same replacement only. Indeed, turning to 
the explicit expression for F given in the Appendix [Xj it would certainly be wrong to use 
there the replacement A — )• A^ (or J(g^) — t- Am(g^)/g^). Instead, the correct procedure is 
that outlined above: the vertex F^ is indeed obtained from the expressions in the Appendix, 
by replacing J — t- J^ (but with no explicit mass terms); all explicit mass terms are next 
added through the totally longitudinally coupled non-perturbative vertex V . 

Actually, it is interesting to ponder about what would happen if one were to introduce 
the gluon mass through the (wrong) procedure of identifying the vertex F' by the simple 
replacement A — > A^ carried out inside F. In such a case one would conclude (after 
some steps) that the self-consistency of the theory would force m?{q^) to vanish identically. 
The precise way how this "self-correction" takes place is intimately related to the so-called 
"seagull identity" [l5| , and will be discussed at the end of the Section |Vl 

IV. GENERAL FEATURES OF THE GLUON MASS EQUATION 

Let us now consider the gluon SDE of Eq. (12. 61) under the light of the analysis presented 
in the previous section. After dynamical gluon mass generation has taken place, one needs 
to consider the modified SDE, which is obtained from (12.61) after (z) replacing the A~^ 
appearing on the Ihs with the A~^ of Eq. (13. ip . and (zz) replacing A — )> A^ and F — )• F' 
inside the integrals of the rhs (see also Fig. [5]). 
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FIG. 5: Diagrammatic representation of the gluon one-loop dressed diagrams before and after 
dynamical gluon mass generation has taken place: the propagators and vertices on the rhs have 
now become massive. 

From this new SDK one can obtain two separate equations, the first one governing the be- 
havior of Jmio^) [to be later involved in the definition of the effective charge, see Eq. (16.171) 
and the second one describing the dynamical mass m^(g^). The general idea is the following: 
the terms appearing on the rhs of the SDE may be separated systematically into two con- 
tributions, one that vanishes as g — )■ and one that does not; the latter contribution must 
be set equal to the corresponding non- vanishing term on the Ihs, namely — m^(g), while the 
former will be set equal to the vanishing term of the Ihs, namely q^Jm{<f)i the so-called 
"kinetic term" . 

Specifically (taking the trace of both sides of (12 .Gp to eliminate the Lorentz indices), the 
rhs may be schematically cast in the form 

gV™(g2) - m\q^) = q^ [l + lCi{q\ m\ Am)] + K,2{q\ m\ A^), (4.1) 

such that g^/Ci(g^, m^. Am) — t- 0, as q^ — > 0, whereas /C2(g^, m^, A^) 7^ in the same limit. 
Thus, for example, a term of the form q^ j^/S.mik)/S.m{k + q) contributes to /Ci, whereas 
a term of the form m'^{q^) j^Am{k)Arn{k + q) should be assigned to /C2- Then, the two 
equations determining Jm{q^) and m?{q^) will read (still Minkowski space) 

Jyniq^) = 1 + /Ci(g^m^ Am), 

m\q^) = -/C2(g^m^Am). (4.2) 

Of course, there is an obvious subtlety that must be addressed at this point. Specifi- 
cally, one may easily envisage the possibility of a term that approaches zero as (g^)", with 
< a < 1. In this case, given that we must factor out a g^ in order to obtain the equation 
for Jm{<f)i such a term would furnish an IR divergent contribution to Jmi.'f)- This would 
be an undesirable feature, given that the Jmil'^) is intimately related to the effective charge 
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of QCD, which is beheved to be finite. The way to treat such a possibihty is to state that, 
should such a term appear, it ought to be directly allotted (in its entirety, without factoring 
out a g^) to the equation for m?{q^). The presence of such a term in the mass equation will 
not affect the value of the mass at q^ = 0, but will in general affect the shape of the result- 
ing curve. Keeping this mathematical possibility in mind, let us point out that the terms 
emerging in the analysis of Section |V] have a very characteristic structure [see Eq. fl4.15p 
and the related discussion], and, at least for them, the scenario contemplated above [ (g^)", 
with < a < 1] is not realized. 

There is an additional point related to the mass equation, which is instrumental for the 
self-consistency of the entire approach. Specifically, a crucial condition for the mechanism 
of dynamical gluon mass generation, developed in a series of articles 
cancellation of all seagull- type of di verg ences, i.e., divergences produced by integrals of the 
type J^Alk"^), or variations thereof 15|. The precise cancellation of such terms proceeds by 



10-13 



16| . is the 



means of the identity [li 

leA'^ik) + ^lAUk) = 0, (4.3) 

where the "prime" denotes differentiation with respect to A;^, i.e., A'^{k) = dAm{k'^)/dk'^. 
Thus, all the ingredients entering into the SDE (most importantly, the vertex) must be 
such that, after taking the limit of the SDE as g^ — ?■ 0, all seagull-type contributions must 
conspire to appear exactly in the combination given on the Ihs of Eq. (14. 3p . 

The relevance and function of the identity (14. 3 p becomes evident when we consider the 
term I{q), given by 

- Uiq) = jeA^{k)A^{k + ^) (fe + g)y»^(_fe+^g^-^fe^^»^(fc) +cj^A^{k), {A.A) 

with c (for the moment) an arbitrary real number. This term appears naturally in the PT- 
BFM framework, and in fact we will find it in the case of the Landau gauge studied in the 
next section. 

Using Eq. (13. ip one may then rewrite I{q) as 

I{q) = h{q) + h{q)^ (4.5) 

with 

2A„(A; + g) - A„(A;) 



-"■fa)--/^^ -'l,-^.;" -^]^^^ 
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c+^j y^A^(fc), (4.6) 



and 



- iiM = l^k^AUk)A,„(k + ,)!!!!|±|j^. 



(4.7) 



In order to establish how the above terms must be assigned among the /Ci and /C2 introduced 
above, let us now take their limit as q^ — )■ 0. Carrying out the appropriate Taylor expansions 



[see Eq. iKWf 



one finds 
-^/i(0) = 



feA'^ik) + ^ fAUk) 

where in the second step we have employed Eq. (14.31) . and 



d 

^+2 



Arr^k) 



(4.8) 



-z/.(0)= feAUk)[m'{k)r. (4.9) 

Jk 

Thus, according to the rules introduced above, the contribution of I{q) to the kinetic term 



IS 



iht{q) 



A;' 



,Am{k + q) - Am{k) d 



{k + qf-k^ 



+ - Am{k), 



(4.10) 



given that /kt(0) = 0, while the contribution of /(g) to the mass equation is 

m?{k + q) — 'm^{k) 



-il^2{q)= I k'^Am{k)Am{k + q)'- ., , .3 12 
Ik [k + qr - k^ 



c + 



f)l^-(^)- (4-11) 



It is clear now that the second term on the rhs of (14. lip is quadratically divergent (and of 
the seagull type). The only way to avoid this divergence is if the coefficient multiplying 
j^Am{k) vanishes, i.e., if c = —d/2. It turns out that, by virtue of the PT-BFM Feynman 
rules, and the fact that gauge invariance is preserved at every level of this approximation, 
the coefficient c comes out precisely equal to ~dj2\ we emphasize that this result can be 
realized only within the PT-BFM framework. Thus, after the seagull cancellation, one is 
left with the first term only, which is perfectly convergent, provided that the mass decreases 
in the deep ultraviolet. As we will see in the next section, in the Landau gauge this term 
accounts for the bulk of the gluon mass equation. 

Even though the term /kt(g) of Eq. (I4.10p does not appear in the rest of our analysis, it 
is important to gain some further intuition on its structure and its behavior for small values 
of g^, especially in the light of the discussion following Eq. (14. 2p . 
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To this end, let us introduce spherical coordinates through the definitions q^ = x, 
k"^ = y, {k + q)^ = z; we then have that z = y + x + 2y/xycos9, and we define 
w = {k + qY — k"^ = z — y = X + 2y/xycos9. The d- dimensional integral measure will read 
in this case 



1 7r2 /,„..n„/ , d 



„ {2nrr{'-^)jo JO 

and we finally recall the elementary integral 



d^sin''-'^ / dyy^-\ (4.12) 



r(ia±i)r(ii±i) 

■ m+n+2'^ 1 '^ even 



d^sin'"^cos"0= <^ n^'^r^) (4.13) 

[ 0, n odd 

It turns out that the /kt(Q') of Eq. (I4.10p may be expanded systematically as a power 
series in q^. To see this in detail, we consider the Taylor expansion of an arbitrary finite 
function f{z) around w = 0, given by 

^^^^^^ = nv) + ^f'iy) + f rfe) + ... (4.14) 

w 2! 61 

where the primes denote differentiations with respect to y (evidently we are assuming finite 
derivatives in the origin). Then, under the integral sign on the rhs of Eq. (I4.10p one must 
collect pieces of a given order in g^ from the various powers of w, using (I4.13p . 

It is clear that when the term f'{y) on the rhs of (I4.14p is inserted into the integral, it 
generates the seagull identity (14. 3p : all the remaining terms will be proportional to positive 
powers of w, and thus, /kt(0) = 0. For example, the g^ term in this expansion is obtained by 
appropriately combining contributions proportional to f"{y) and f"'{y)- Using again (I4.13p . 
after a sequence of partial integrations, we find 

lUq) = ^ (^^) J^ eAlik) + 0{q% (4.15) 

In order to check the validity of Eq. (I4.15P let us compute I\<±{q) for the simple case of a 
massive propagator with a "hard" (momentum-independent) mass 

A^(g) = -Y^- (4.16) 

The integrand in the first integral on the rhs of Eq. (14.1 ip simplifies to 

^ A^{k + q)-Am{k) fc2 

{k + qf-k^ {k^ - m^)[{k + qY - m^X ^ ' ' 
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Then, using the dimensional regularization identity 



2m2 



i(raF-(^-^)ira- ("«) 



w^) = i6^/'^-'"r + ^i^i- ^'-''^ 



it is relatively straightforward to demonstrate that 

m^ /"^ 1 1 f 1 : q^x{x — 1 

'0 

Evidently, /kt(0) = 0, as expected, while the expansion of the logarithm furnishes immedi- 
ately the result 

4t(g) = -^^+0(g^). (4.20) 

On the other hand, substitution into the general formula (I4.15P of the propagator in 
fliA6|l yields 

In dimensional regularization, around rf = 4, we have that d = 4 — e, and therefore, only the 
divergent part of the integral contributes, i.e., 



U,) = -4 (^ ) 2 



id\ (2 



167r2 V 4 / V e 



0{q' 



1 q 



2 



+ 0(g^), (4.22) 



IGtt^ 6 

which indeed coincides with (14.201) . 

Notice finally that the main contribution to the kinetic term does not originate from 
/kt(q'), but rather from a term of the form 



2 ,"Am(fc + g) - Am(/c) 



^1 ";..;i^-p ■ <^-^^' 

which, for the simple massive propagators of fl4.16p may be easily calculated, giving rise to 
the standard logarithmic correction associated with the RG, with the additional feature of 
being IR safe due to the presence of the mass in the argument of the logarithm. 

V. THE GLUON MASS EQUATION IN THE LANDAU GAUGE 

We now proceed to the actual derivation of the explicit form of the mass equation in the 
Landau gauge. Specifically, in this gauge the rhs of Eq. (12. 6p will be given by the terms Ai 
listed in Eq. (I2.16p . where now we must carry out the replacements A — )• A^ and F — )► F'. 
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Following the rules explained in the previous section, and defining 

A,(g) = Tr[Ar(g)], (5.1) 

the mass equation is given by 

"^^^f- 'd-1 [1 + G(g2)]2 ' ^^■^> 

and therefore, one should determine the closed form of the quantities [y4j(g)]„2. 

There is a simple observation, particular to this gauge, which simplifies the entire pro- 
cedure considerably. Specifically, in the Landau gauge, the derivation of the gluon mass 
equation does not require the knowledge of the closed form of the vertex V, which captures 
the effects of the massless bound-state poles. 

To see why this is so, let us first note that the vertex V appears only in the terms A'^''{q) 
and A'^'^{q), the only place where the replacement F — t- F' may be carried out. Given that 
the vertex T^^^^ appearing in the term A'^'^{q) is contracted by {k + q)'^, the result of this 
operation is the STl satisfied by F', namely 

p^T' = Fip') [A„i(r2)P;(r)#,„(r,p,g) - A;^\q')P^iq)H,,iq,p,r)] . (5.3) 



whose validity assumes the existence of V but does not depend on the details of its closed 
form. 

As for the term A^'^{q) one starts by noticing that (i) the V is already contracted by 
two projection operators P°'P{k)P'^"{k + q) and (ii) since in the PT-BFM formulation the 
truncated 11^'^ (g) (defined in terms of Ai — A^) is transverse, one may contract both sides 
of Eq. fl2.15p by the projection operator P^^i^q) for free, i.e., write 

5 

= /CA$^Ar'(g)P;(g). (5.4) 

The main effect of this operation, as far as the term A'l^{q) is concerned, is to trigger 
Eq. fl3.4p . and so, all explicit reference to V vanishes. 

In order to forestall any possible confusion, we hasten to emphasize that one should not 
conclude from the above argument that the existence of the vertex V is irrelevant for the 
entire construction. On the contrary, the vertex V is crucial for the implementation of this 
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particular approach. In particular, if the V did not exist {i.e., if it were vanishing identically) 
the WI of (13.51) would be invalidated, and, as a result, Il^'^{q) would fail to be transverse, in 
which case, evidently, one could no longer contract both sides of Eq. (12.151) by the projection 
operator Pl^/{q) for free. 

We can now proceed with the actual calculation. It is clear that the term A'^'^{q) cannot 
possibly contribute to the mass equation, since A'^'^{0) = 0. Furthermore, with the exception 
of A'j^'^{q), which will yield a direct contribution, the remaining three terms A'^'^{q), Ai^'^lq), 
and A'^'^{q) contribute to the mass equations an amount that arises as the deviation from 
the seagull cancellation, i.e., they furnish a term analogous to the I2 of Eq. (14. 7p . 

To see this, let us first retain the contributions of these three terms that survive individ- 
ually the g^ — )■ limit. The term Ai reads 

1 



Aiiq) 



^^Lp:{k)p.'{k + q)T'::''^PMAUk)AUk + g). 



(5.5) 



Now, using for the vertex Tm the tensor decomposition (1A2|) with (a, /i, u) -^ (z/', p, a) and 
r ^ k, p ^ —k — q, it is straightforward to establish that the tensorial structures £2, £5 and 
is will be annihilated by the transverse projectors appearing in (15. 5p . while, ignoring terms 
that will again vanish due to the transverse projectors, ii, £3, £7, and £g are at least of order 
q. Finally, since £10 = 0, we find the result 



T^i^'^'^iq, k, -k-q) = 2k^'r [X,{q, k, -k - q) + k^X,{q, k, -k - q)] + 0{q) 

.{k + qfJm{k + q)-eJm{k) 



2r c/^"- 



{k + qf - F 



In addition, since 



^(0) 



TZ3{q.k,-k-q)k^Pi:.{q) = 2g, 



aP 



e- 



{k-qf 



0{q). 



+ 0{q). 



(5.6) 



(5.7) 



we finally obtain 
A^{q) = 2{d-l) 



{k-qf 



{k + qfJm{k + q)- k^Jm{k) 



{k + qf - A;2 



A„(A;)A„(fc + g) + 0(g). 

(5.8) 



Similarly, from A2^{q) we obtain 

A2{q) = - 
while ^4^(5) contributes simply 

A,{q) 



d-2 + 



{k-qf 



{d-lf 



k ^m(^j 

(A; + g)2 ■ 



d 



Am{k). 



(5.9) 



(5.10) 
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The terms in Eqs. (15. Sp . (15. 9p and (I5.10p are individually non- vanishing as q^ — )■ 0, but 
their final contribution to the mass equation is controlled by the seagull identity, which forces 
a large part of their sum to vanish, thus reassigning them to the kinetic term. Specifically, 
if we use Eq. (13.11) to substitute the terms containing J^ in the numerator of the integral 
on the rhs of Eq. (15.81) . [i.e., k'^J{k) = A^{k) + m^(A;)] the sum of these three terms gives 



[Ai + A2 + Ai]{q) = [Ai + A2 + A4]kt(g) + [^1 + A2 + Ai]m2{q), 



(5.11) 



where 



[A + A2 + ^4]kt(g) = -2{d-i] 



d-2 + 



and leaves as residual contribution 



2 {k ■ qf] A^{k + q) - A^{k) 



F- 



{k-qy 



{k + qy-k^ 



k'AUk) (d-l) 



k'^q^ J {k + qy 



d 



A^{k), (5.12) 



[Ai+A2+A^)]^2{q) = 2{d-l) 



2 {k ■ qy] m'^{k + q) — m^{k) 



k'- 



q' 



{k + qy - P 



A^{k)A^{k+q). (5.13) 



It is now easy to verify that, by virtue of the seagull identity, the rhs of (I5.12p vanishes 
as g^ — 7- 0. Indeed 



[A + A2 + A4]kt(0) = -2(rf-l) [sinHeA'^ik)- I 

Jk Jk 



and using that [see also Eqs. (14. 12p and (14.130 above] 



(d-l)- sin^ e 



jd-l] 
d 



sin' ef{k) = ^ [f{k), 

k 0- Jk 



Am{k), 
(5.14) 

(5.15) 



it is elementary to demonstrate that the rhs is exactly proportional to the expression on the 
rhs of Eq. (14.31) . and therefore vanishes. 

We next consider the term ^3. After taking the trace we find 



A,{q)= P^^{k) 



(* + *)'r,^P:'fe,^„„). 



(5.16) 



Ik (k + qy 

When inserted into the expression for A^i^q), the first term on the rhs of (13. 9p will give the 
result 

r Flh -U n\ 

(5.17) 



^^^ + f a(g,-fc-g,fc) + 0(g). 



(A; + qy 
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which contributes to the kinetic term, since in the g — )■ hmit vanishes due to the second 
identity in (IA5|) . which in this hmit gives J36| 

a{Q,-k,k) = F-\k). (5.18) 

The second term on the rhs of (13.91) yields instead a surprisingly simple contribution to the 
mass equation. Specifically, using the definition fl2.1ip we obtain 



= m^(g^)%^P'^^(g) 

= i^rh'{q')Giq'). (5.19) 

g Ca 

On the other hand, the second of the background quantum identities (I2.14p implies (see also 
Appendix [B] for an alternative derivation of this result) 

ffi\q^) = [l + Giq^)]m\q^), (5.20) 

so that one finally finds the contribution 

[AsUiq) = t^Giq')[l + Giq')]m\q'). (5.21) 

The next step is to substitute the above results on the rhs of the mass equation of 
Eq. (15. 2p . In doing so, we move to the Euclidean space, by setting J^= ij^ and q^ = — g^, 
and using 

AM) = -M-qD; rnliql) = m'i-ql); G^l) = G{-ql). (5.22) 

Then, from Eq. (I5.13P and (I5.2ip . we arrive at the final form of the mass equation, namely 



2/ 2\ 2(7 Ca 



A;^- 



{k-q) 
g2 
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m?{k + q) — m^{k) 



A^{k)AUk + q), (5.23) 



(fc + g)2-A;2 

where we have suppressed the suffix "E" . 

Finally, we are now in position to address the question posed at the end of Section IIIIl 
namely what would happen if we were to introduce the gluon mass by the simple replacement 
g^J(g^) — 7- A,^^(g2) carried out inside F, i.e., without resorting explicitly to the vertex V 
(with the crucial properties assigned to it). The basic observation is that the main bulk of 
the mass equation, namely the rhs of Eq. (I5.13p . emerges as a residual contribution that 

22 



survives the seagull cancellation. However, within this hypothetical scenario, the term ^1(13') 
in Eq. (15.81) would be instead given by 



A(g) = 2(rf-l 

= -2(rf-l) / 

Jk 



^,_,.. ^,^]^-^\k + q)-/\-^\k) 



.2 (fe • ^f 
-,2 



I 

k 



2 {k-qy 



q2 



{k + qf-k^ 
Am{k + q) - Am{k) 



Am{k)A^{k + q) 

(5.24) 



{k + qy-k^ 

thus, participating in the cancelation of Eq. (I5.14p . as before, but leaving no residual con- 
tribution, i.e., [Ai + A2 + A^lm^lq) = 0. Then, the only contribution to the rhs of the mass 
equation would be that of [A3]m2(g) in Eq. (I5.2ip : this contribution would be still there, 
because within this alternative scenario the full vertex F' is still assumed to satisfy the full 
STIs of Eq. (13. 9 p [but with no reference to V]. Therefore, the resulting mass equation [the 
equivalent of Eq. (I5.23p ] would read 

which would simply imply m?{q^) = 0, i.e., no dynamical mass generation. 

VI. NUMERICAL ANALYSIS 

In this section, we will first derive an approximate version of the mass equation (I5.23p . 
which will facilitate the numerical treatment while retaining the main features of the full 
equation. Then, using as input for the functions A(g^) and F{q^) [appearing in (I5.23P ] the 
available lattice data, we solve the equation numerically for the gauge groups SU{2) and 
SU{3), thus obtaining the (approximate) form of im?{q^). Then, using Eq. (13. ip . together 
with the A(g^) of the lattice and the m?'{q^) obtained from the mass equation, we will 
extract the (approximate) form of Jm{q^). As a basic application, these ingredients will be 
subsequently combined to form the gluon mass entering in the RG-invariant combination 
associated with the definition of a non-Abelian effective charge. 

A. Approximate version of the mass equation 

We now proceed to the analysis of the mass equation (15.231) . The difficulty in dealing 
with this equation in its full version resides in the fact that the unknown function m? 
appearing on the rhs depends on both the angular and the radial coordinates {6 and y, 
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respectively). To circumvent this problem we will employ certain standard approximations, 
in order to eliminate the angular integration. However, before embarking into the derivation 
of the approximate version of (15.231) , we can extract useful of information about the global 
behavior of m^ from its g^ — )■ limit. 

Specifically, let us employ the notation introduced in ( I4.12p . and consider the limit of 
Eq. f l5.23p as g^ — )■ 0. Since it is known that L(0) = in four dimensions [41], Eq. f l2.13p 
implies that 1 + G'(O) = F"^(0), so that we get 

"^'(0) = Ig'CAFiO) [e[m'{k)]'A\k) 
2 Jk 

= -3/CaF(0) [m\k)A{k)[eA{k)]' . (6.1) 

Jk 

Obviously, in the kernel of above equation there is no dependence on 6, so that the angular 
integral can be done exactly, and one is left with the final equation 

m'iO) = -^asFiO) j dym'iy)[y'A'iy)]\ (6.2) 

where as = g'^/4:7c and, as usual, y = k"^ (the prime indicates now derivatives with respect 
toy). 

Equations (16. ip and (16. 2p furnish a rather interesting constraint on the structure of the 
full gluon propagator. Indeed, it is clear that due to the positive sign in front of the first 
line of Eq. (16. ip . solutions of (I5.23P leading to a positive m'^{0) cannot be monotonically 
decreasing; or, seeing it from the point of view of Eq. (16. 2p . the kernel [y'^A'^{y)]' must 
reverse sign and display a "sufficiently deep" negative region at intermediate momenta, in 
order to obtain m'^{0) > 0. This is a highly non-trivial requirement, because, to the best 
of our knowledge, there is no a priori fundamental reason why the full gluon propagator 
propagator should show this particular behavior. 

We now proceed to the derivation of an approximate version of (I5.23P that will reproduce 
in the g^ — )■ limit Eq. (16. ip . and therefore implement the important constraint that this 
equation entails. 

Let us then denote by R{q) the integral appearing on the rhs of (I5.23p : using the simple 
identity 

{k-qr = \{[{k + qr-kr-2q'[{k + qr-e] + {qr}, (6.3) 

we see that the second term above, when inserted back into R{q), vanishes upon integration, 
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and therefore one is left with 

R{q) = Ri{q) + R2{q), (6.4) 

where 

2 q^^ m?{k + q) —ni?{k) 



R2{q) = -7^ jm\k)[{k + qf - A;2]A„(A;)A„(A; + q). (6.5) 

^q Jk 

To cast Ri{q) and R2{q) into a form suitable for solving the corresponding dynamical 
equation, we first introduce the by now familiar spherical coordinates and then split the 
radial integration into two intervals 



OO /"X fOO 



dy = dy+ dy, (6.6) 

Jo Jx 

SO that in the second integral since y > x always, we can expand the integrand according 
to fl4.14p . Proceeding in this way, and observing that partial integration gives 



dyy'[m\y)rA\y) = -m'{x)x'A\x) - / dym\y)[y'A'{y)]' (6.7) 

^ X J X 

we obtain 

167T^R,{x) ^ A(x) Pdyy (y - -] ""'^"^^ ~ ""'^^^ A(y) - m\x)x^A\x) 
Jo ^ 4/ x-y 

- / dym'iy)[y'A'iy)]\ 

J X 

IQ-K^R^ix) ^ ]^j dyym\y)(l-^A\y) + -^j dy m\y)[y^ A^y)] . 



(6.8) 



Finally, since as shown in 4l|, |42|, L{x) is considerably smaller than G{x) in the entire 
range of (Euclidean) momenta, we can use the approximation 1 + G{x) ~ F^^(x); thus, we 
obtain the approximate equation 



2/ A 2fr.^F{^) . OisCA 



m\x) = m\0)j^ + -^ F{x)R{x), (6.9) 



with 



R{x) = i£dyym^(y)(l-|)A^(y) + A(x)^%y(y-|) ^'^^]_^'^^^ A(y) 

- m\x)x^A\x) + ^-j dym\y)[y^A\y)~^, (6.10) 

and m^(0) given in Eq. (16. 2p . Evidently, i?(0) = 0. 
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Gluon propagator in SU(3) 

=5.7, ii = 4.3GeV 
c V=64' 
□ V=72* 
A V=80' 
fit 




q[C3eV'] 



Gluon propagator in SU(2) 

p=2.2, II = 2.2 GeV 

o V=128' 
fit 




q [GeV] 



FIG. 6: Lattice results for the SU{3) (left) and SU{2) (right) gluon propagator, renormalized at 
fj, = 4.3 GeV and fi = 2.2 GeV respectively. The continuous lines represents our best fits to the 
data obtained from Eq. (j6.1ip . 



B. Lattice ingredients: Gluon propagator and ghost dressing function 

The two main ingredients of the mass equation (16.91) are the gluon propagator A(g^) 
and the ghost dressing function F{q^). Of course, A(g^) is composed by J{q^) and m?{q^), 
as dictated by Eq. (13. ip . but, as mentioned in the Introduction, the derivation of the cor- 
responding equation for J{q^) is beyond our powers at this point, mainly due to lack of 
.cnowledge of certain of its ingredients. Similarly, F{q^) satisfies its own SDE (see, e.g., 
4l|). which would furnish yet another equation in a complicated coupled system. For the 
purposes of the present work, which is the preliminary scrutiny of the mass equation (16. 9p 
appearing for the first time in the literature, we will instead resort to the high quality lattice 
data available, and use them as inputs inside (16. 9p . 

In order to do that, we start by showing on the left panel of Fig. [6] the lattice data for 
A(g^) obtained in 6|, corresponding to a S'f/(3) quenched lattice simulation, renormalized 
at yU = 4.3 GeV; on the right panel of the same figure, we show the quenched SU{2) lattice 
data obtained in [l|, renormalized at /i = 2.2 GeV. 



As has been discussed in detail in the literature Uil, \Va 



16|, both sets of 



attice data can 



be accurately fitted in terms of a IR finite gluon propagator of the form 43 1 

^-\q') = M\q^) + q' 



\^l^CAgi^^[q' + PiM\q^^ 



967r2 



^^- 



(6.11) 
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Ghost dressing function SU(3) 

|J=5.7, n = 4.3GeV 
o V=64* 

A V=' 
fit 




Gliost dressing function SU(2) 

p=2.2, n = 2-2 GeV 

o V=128' 
fit 




q[C3eV'] 



q [GeV] 



FIG. 7: Lattice results for the SU{3) (left) and SU{2) (right) ghost dressing function, renormalized 
at ^ = 4.3 GeV and /u = 2.2 GeV respectively. The continuous lines represent our best fits to the 
data obtained from Eq. ()6.13p . 



where 49 1 



2/'^2n 



M\q 



rrin 



g2 + p^ml 



(6.12) 



The function M'^{q^) controls the value of A~^(g^) at the origin; evidently, 
A~^(0) = M^(0) = m^l p2. The best fits (shown by the continuous lines in Fig. [6]) corre- 
spond to the following values of the fitting parameters: 

• 5t/(3) case: m^ = 520 MeV, gf = 5.68, pi = 8.55, ps = 1-91; 

• SU{2) case: rriQ = 867 MeV, gf = 10.80, pi = 1.96, pa = 2.68. 

Turning next to the ghost dressing function, on the left panel of Fig. U\ we show the 
SU{3) lattice results of [6[, renormalized as before at p = 4.3 GeV; on the right panel we 
plot instead the results for the SU{2) case 

clearly seen, both functions saturate in the deep IR at the constant value 
can therefore be fitted in terms of the expression 



l[, renormalized at p = 2.2 GeV. As can be 



12 



50 



\q' 



1 I ^Cp^gl^^fq^ + p,M\q'^ 



4 487r2 V p2 



5l|, and 
(6.13) 



with M'^{q^) given by Eq. (I6.12p . but changing the parameter p2 — ?■ P4. 
The best values for the fitting parameters are: 

• 5f/(3) case: gl = 8.57, m = 520 MeV, pa = 0.25, p^ = 0.68; 

• SU{2) case: g^ = 15.03, m = 523 MeV, pg = 0.21, p^ = 0.78. 
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q [GeV] 



q [GeV] 



FIG. 8: The kernel [q^ I!^'^ {q^)]' appearing in Eq. ^^ obtained from the 5[/(3) (left) and SU{2) 
(right) lattice data. In both cases one clearly sees the behavior expected for getting a positive 
value for m?{Q). The zero crossing happens at q^ ^ 0.85 and q^ ^ 1.1 respectively. 

C. Solutions of the mass equation and extraction of Jm{q^) 

After presenting the precise form of A(g^) and F{q'^), the next task is to find solutions 
of the approximate mass equation (16. 9p . 

To begin with, we compute (for both gauge groups considered) the derivative of the 
gluon dressing squared, [y^A^(t/)]', entering into the condition (16. 2p . As mentioned earlier, 
the behavior of this quantity provides a rather direct criterion for the existence or not of 
positive-definite mass solutions, and in particular m'^{0) > 0. Specifically, the absence of a 
negative region from this derivative immediately excludes such solutions, while a relatively 
shallow "well" makes their existence unlikely. 

In the results shown in Fig. |8] we clearly see that both derivatives change their sign 
in the intermediate momenta region, which, as previously explained, constitutes precisely 
the required behavior. This behavior is to be contrasted with that of simple propagators, 
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which fail 



such as l/iq"^ + m"^), or the Gribov-Zwanziger propagator q'^/{q^ + m^ 
to provide the necessary negative region (in fact the derivative is positive everywhere). It 
should be noted that, instead, the "refined" version of the Gribov-Zwanziger propagator [52 1 
is expected to furnish a considerable negative region, given that it is known to provide a 
good fit to the lattice data. 

Of course, the aforementioned criterion can only serve as a necessary but not sufficient 
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Gluon mass in the Landau gauge - SU(3) 

a (n) = 0.59 and m'(0)=0.14 
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Gluon mass in the Landau gauge - SU(2) 

aju) = 3.2 and m'(0)=0.28 
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FIG. 9: The solution for m?{q^) obtained through the approximate mass equation ()6.9p for SU{'i) 
(left) and SU{2) (right). 

condition: to get a positive definite value for m^(0) one still needs to demonstrate that the 
negative region q^ > q^ (with Qq the value where the curve is zero) furnishes more support 
to the integral of Eq. 06. 2p than its positive region. 

To proceed with the actual determination of m^{x) from Eq. (16. 9p . we substitute the 
quantities A{y), F{y) and Ca for the SU{3) and SU{2) gauge groups and solve for the 
unknown function. In both cases the value of m^(0) is a boundary condition, fixed through 
the value of the corresponding lattice gluon propagator at the origin, i.e., ?7i^(0) = A^^(O). 
Specifically, for 5f/(3) we have that A-i(O) ^ 0.14 while for SU{2) A-i(O) ^ 0.28 

The solutions obtained are shown in Fig. [HI the values for a^ needed to satisfy the 
boundary condition are as = 0.59 and as = 3.2 for SU{3) and SU{2) respectively. Notice 
that the masses corresponding to both gauge groups display the same qualitative behavior, 
and, as expected, are clearly non-monotonic functions of the momentum. 

From the solutions for m^(g^) obtained above, and the lattice results for A(g^), we may 
now extract the approximate form of the "kinetic term", Jm(g^). Specifically, JmiQ"^) can be 
determined (in Euclidean space) through Eq. (13.11] , namely 



Jm{q 



(6.14) 



Notice that special care must be taken in the g^ — > limit of Eq. (I6.14p . In the region 
of small momenta, Eq. (I6.14p has a delicate cancellation between the denominator and the 
numerator, which also tends to zero in this limit, since A~^(g^) — > m'^{0). In order to avoid 
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SU(2) 

extrapolation 
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FIG. 10: Values of Jmiff') obtained from Eq. ()6.14p (white circles) using the SU{'i) gluon propa- 
gator and the corresponding extrapolation towards the g^ — t- limit (continuous line). As usual 
we show both the 5C/(3) (left) and the SU{2) cases. 

spurious distortion in the IR behavior of Jmi^f)-, we will extract Jmiff) until certain (small) 
value of q^ past which we will do an extrapolation towards q^ — )■ 0. The results of this 
procedure are shown in Fig. [TOl where, for both the SU{?>) (left) and the SU{2) (right) 
cases, we display the points obtained directly from Eq. (I6.14p as well as our extrapolation 
curves. 

Knowledge of m?{q^) and Jm{q^) allows one to determine the approximate form of the 



(formally) RG-invariant 



^ 



uon mass that appears naturally in the definition of the QCD 



effective charge [l6|, l2l|, |22|, |53|. Let us recall that, due to the Abelian Wis satisfied by the 
PT-BFM Green's functions, the propagator A(g^) absorbs all the RG logarithms, exactly 
as happens in QED with the photon self-energy. As a result, the product 

do(g^) = glMq') = /A(g2) ^J (g^), (6.15) 

forms a RG-invariant (/^-independent) quantity. As has been explained in the recent litera- 
ture 15|, d (g^) may be cast in the form 



d{q' 



q^ + m?'{q'^) 



(6.16) 



with 



m\q') = m\q')J-\q'). 



(6.17) 
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FIG. 11: The RG-invariant mass, m'^{q'^), defined in Eq. dOOjl for tlie SU{3) (left) and SU{2) 
(right) cases. 

The two factors defined above are individually RG-invariant; the dimensionful quantity 
corresponds to a massive propagator with a momentum dependent mass, while the dimen- 
sionless factor ^^(g^)/47r defines the effective charge. 

Next, using the BQIs (12.141) to relate the components of A(g^) to the corresponding ones 
of A(g^), we get 



JUq') = [l + G(g2)]V„(g2), 
m\q') = [l + G(g2)]W(g2), 



and therefore 



which finally furnishes the relation 



-l^„2\ 



m'{q') = m\q')J;;,\q 



(6.18) 



(6.19) 



(6.20) 



We are now in the position to determine the mass rff'{q^) by simply forming the ratio 
of the plots presented in Fig. [9] and [TOl The result is shown in Fig. [TTl as can be seen, 
in the S'f/(3) case m^(g^) corresponds roughly to a monotonically decreasing function (see 
also 541), with m(0) ^ 580 MeV. Finally, for the SU{2) case we obtain m(0) ^ 480 MeV. 
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VII. DISCUSSION AND CONCLUSIONS 

In the present work we have derived the dynamical equation that determines the evolution 
of the gluon mass in the Landau gauge, using as our starting point the "one-loop dressed" 
SDE for the gluon propagator in the PT-BFM scheme. The entire construction hinges on the 
crucial assumption that a special vertex, denoted by V, is dynamically generated, according 
to the philosophy and formalism associated with the Schwinger mechanism. The role of 
this vertex is to maintain gauge invariance (as expressed through the STIs satisfied by the 
Green's functions of the theory) in the presence of a dynamical mass. Interestingly enough, 
the derivation of the mass equation does not depend on the specific closed form of that 
vertex. 

The equation for the gluon mass derived here, given in (15.231) . and in particular its limit 
in the deep IR, imposes a rather strong constraint on the form of the full gluon propagator in 
the region of intermediate momenta of about (1-5) GeV^. In this specific range of momenta 
the shape of the gluon propagator must be such that the derivative of the square of the 
gluon dressing function [g^A^(g^)]' becomes sufficiently negative, thus ensuring eventually 
the positivity of the gluon mass. 

We emphasize that the central result of this article, Eq. (15.231) . does not exhaust all pos- 
sible contributions to the gluon mass equation. Specifically, Eq. (I5.23P captures only the 
part of the equation originating from the "one-loop dressed" gluon SDE. In order to deter- 
mine the corresponding contribution coming from the "two-loop dressed" gluon SDE one 
must identify the seagull cancellation mechanism (and the corresponding "seagull-identity" ) 
that operates at the "two-loop dressed" level. The identification of the "two-loop dressed" 
analogue of Eq. (14. 3 P requires (among other things) some very specific information on the 
structure of the four-gluon vertex, at least in the special kinematic limit of vanishing external 
momentum. Calculations in this direction are already in progress. 

It is important to warn the reader about some additional limitations afflicting the present 
work, related to the renormalization properties of the mass equation, and the dependence of 
the various quantities, most importantly of the gluon mass, on the renormalization point /i. 
When dealing with the mass equation of Eq. (15.230 we have tacitly assumed that the mul- 
tiplicative renormalization has been carried out, thus rendering all quantities finite (but fi- 
dependent). In carrying out the SDE renormalization one usually resorts to the momentum- 
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subtraction (MOM) scheme; in our case this choice is further motivated by the additional 
fact that this is the scheme employed for the renormalization of the lattice data that are 
used as input into Eq. (15.23 p . Of course, given the gluon mass in the Landau gauge is not a 
RG-invariant quantity, there is a residual dependence on /i, which, in principle, should cancel 
out against analogous contributions when a RG-invariant combination is formed (this type 
of powerful cancellation has been presented in J4l| for the QCD effective charge). However, 
the approximations employed in the process of the renormalization may distort the exact 
dependence on /i. Specifically, the renormalized version of Eq. (I5.23P displays a dependence 
on some of the renormalization constants Z involved, as happens typically in the treatment 
of SDEs. This fact in itself is normal, but makes the further treatment ambiguous, because 
the correct cancellation of the residual dependence on the UV cutoff (induced by the pres- 
ence of the Z) requires the knowledge (among other things) of the transverse (automatically 



conserved) part of the full vertex F J36|. Therefore, the next step has been to set Z = 1, 
a fact which, in general, is known to alter the dependence of the solution (in this case of 
re? ) on [i. In fact, the situation appears to be very similar to what happens typically in 
the studies of chiral symmetry breaking through the standard gap equation. In this latter 
context, the various approximations associated with renormalization introduce characteristic 
artifacts; for example, the value of the anomalous dimension of the dynamical quark mass 
is distorted, a problem that is usually compensated by modifying accordingly (by hand) the 
kernel of the gap equation. Needless to say, it would be very important to improve on any 
of the above points, but at present this appears to be technically rather difficult. 

Given that the existence of a non-trivial vertex V is of central importance, it would 
be absolutely essential to establish its existence. This can be done following two distinct 
but complementary approaches. First, one may write down the most general longitudinal 
structure allowed by Lorentz symmetry and then use the WI and STIs that the V is supposed 
to satisfy \e.g.., ( 13. Sp ] to actually determine the form of the various form factors, in the spirit 
of |36| . Second, one may address the dynamical question of whether such a nonperturbative 
vertex may be actually produced by the strongly coupled Yang-Mills theory. In fact, the 
main characteristic of the vertex V ^ which sharply differentiates it from ordinary vertex 
contribution, is that it contains massless bound-state poles. In principle, the dynamical 
formation of such poles must be studied by means of a homogeneous Bethe-Salpeter equation. 



following the methodology developed in |3ll-l33| . We hope to be able to pursue some of these 
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points in the near future. 
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Appendix A: Explicit form of the vertex IT and Wm 



The lon gitu dinal part of the vertex F (and therefore also that of Fm) has been con- 
structed in J36l] by simultaneously solving the Ward and Slavnov- Taylor identities presented 
in Eq. fl2.19p : for the kinematics, see Fig. [2J Specifically, the longitudinal part is written as 

10 

r"^'^(g,r,p) = J2M'l,r,p)er{q,r,p), (Al) 



j=i 



in the standard basis £,• of 



39] 



^a^iu 



5 



i'^'"' = (g - r^g'^l^ 

Qi^" = (r - pYg^" 



.q^gf^u 



if-"" = {q- rflq^'r'^ - (q ■ r)g"''] 
(r — p)"[r'^p^ — (r ■ p)g'^''] 



f/a/iu 



t%^^ = (p — q)^\p°'(f — (p ■ q)g'^ 



(A2) 



and the Aj are given by 



Ai(g,r,p) = -J(g2) ^-p^lj^^^p [r^) + [2a,pg + p^fe^^^ + 2(g ■ r)4pg]i^(p^)} 



+ i'^(^') 



q')b,prF{p') +-J{p')p'b,rpF{r 



X2{q,r,p) = -J{q') {{q^ - r^)b,r,F{r^) + [2ar„ + (r^ - q^)brp, + 2(g ■ r)dr„]F{p^)} 



+ ^J{r') [-2 + p\,rF{p')] + \j{p') (r^ - q') b,r,F{r') 



X3{q,r,p) 



q2 _ J.2 



J{<f) [o-rpq - {q ■ p)drpq] - J{r^) agpr - (r ■ p)d, 



iiqpr 



X,{q,r,p) = ij(g')g' [bprgF{r^) + brpgF{p')] + ^J{r^) [2 - q\prF{p' 
+ ^ J(p2) [2 - q%rpF{r' 
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X,iq,r,p) = -J{q^){p^ - r') [V.^(^') + K^Fip')] + ^Ar') U + {r^ - P^)\,rF {p^) 



4 



Xe{q,r,p) 






X7{q,r,p) = Xi{q,p,r) 

Xsiq^r^p) = -X2{q,p,r) 

Xg{q,r,p) = X^{q,p,r) 

Xioiq,r,p) = ^ { J(g2) [V.^(^') - K,F{p)] + J {r^) F {p^)\,r - J {p") F {r^)K,] ■ (A3) 

The functions aqrp = a{q,r,p), etc are the form factors appearing in the tensorial decompo- 
sition of the ghost-gluon kernels H^^{p, r, q) and H^^{p, r, q), namely 



-^ufiyPj ^) q) QpLV^jq^p f f^qu^qr-p ~r qp,PuCqrp ~r qvPpOiqrp ~r PpPuS-qrp) 

and similarly for H. They satisfy the non-trivial all-order constraints 

F{r'^)[aprq - (r ■ p)hprq + {q ■ p)dprq] = F{p^)[arpg - (r ■ p)brpg + {q ■ r)drpq] 



(A4) 



F{r'^)[aqrp - {q ■ r)bqrp + (g ■ p)dgrp] = 1. 



(A5) 



,2/„2\ 



Appendix B: On the relation bet^veen m {q ) and m {q 

In Section |V] we have assumed that the relation (15.201) between the masses m{q) and 
m{q) holds. This is tantamount to claiming that the BQIs (12.141) hold after dynamical mass 
generation has taken place. 

To further substantiate this claim, let us consider the SDE for the QB gluon self-energy 
n. If we keep dressed the background side of the equation, we can still truncate meaningfully 
the SDE retaining only the one-loop dressed gluon contributions, which now read 






IJflU 



]- g'CA fr^'Uq, k, -k - q)A^P{k)A^^k + q)T,^^{q, k, -k-q), 



2)fiu 



g'c. 



g,^ /a^(A:) - /A,.(fc) 

Jk Jk 



(Bl) 



The projection to the Landau gauge gives rise to three terms only, which coincide with Ai, 
A2 and A^ of Eq. (12.161) . Then writing 



A-i(g^) = q'J{q') - m\q') = q' + tU{q 



(B2) 
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it is relatively straightforward to establish that the rhs of the equation for m is determined 
by the the mass term of Ai only. Specifically, using the result (15.131) one has (Euclidean 
space) 



m\q') = 2g'CA 



Ik 
,2, 



F- 



{k ■ qY] m'^{k + q) — m^{k) 



q^ 



{k + qf-k^ 



Am{k)Arn{k + q) 



= j^[A + A, + A,]^, , (B3) 

where in the last step we have used Eq. (I5.13p . 

Substituting the above result, together with (I5.19p . into Eq. (15. 2p . we find 

^2/2) _ m^(g') ^ m\q^)G{q^) 



[1 + G(g2)]2 [1 + G(g2)]2 
m^(g^) 



1 + G(g2) 
namely the relation of Eq. (I5.20p . 
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